arXiv: 1501.00947vl [astro-ph.SR] 5 Jan 2015 


Astronomy & Astrophysics manuscript no. AA_2014_24333 ©ESO 2015 

January 6, 2015 


Bulge RR Lyrae stars in the VVV tile b201^ 

F. Gran^’^, D. Minniti^’"^, R. K. Saito^, C. Navarrete'’^ I. Dekany^’\ I. McDonald^, R. Contreras Ramos^’^, 

M. Catelan^’^ 

' Instituto de Astrofisica, Pontificia Universidad Catolica de Chile, Vicuna Mackenna 4860, Casilla 306, Santiago, Chile 
e-mail: fgran@astro .puc. cl 

^ Millennium Institute of Astrophysics (MAS), Santiago, Chile 

^ Departamento de Ciencias Fisicas, Universidad Andres Bello, Repiiblica 220, Santiago, Chile 
Vatican Observatory, V-00120 Vatican City State, Italy 

^ Universidade Federal de Sergipe, Departamento de Fisica, Av. Marechal Rondon s/n, 49100-000, Sao Cristovao, SE, Brazil 
^ Jordell Bank Centre for Astrophysics, Alan Turing Building, Manchester, M13 9PL, UK 

Received June 04, 2014; accepted November 24, 2014 


ABSTRACT 

Context. The VISTA Variables in the Via Lactea (VVV) Survey is one of the six ESO public surveys currently ongoing at the VISTA 
telescope on Cerro Paranal, Chile. VVV uses near-IR (ZYJHKs) filters that at present provide photometry to a depth of ~ 17.0 mag 
in up to 36 epochs spanning over four years, and aim at discovering more than lO* variable sources as well as trace the structure of 
the Galactic bulge and part of the southern disk. 

Aims. A variability search was performed to find RR Lyrae variable stars. The low stellar density of the VVV tile b201, which is 
centered at {{,b) ~ (-9°,-9°), makes it suitable to search for variable stars. Previous studies have identified some RR Lyrae stars 
using optical bands that served to test our search procedure. The main goal is to measure the reddening, interstellar extinction, and 
distances of the RR Lyrae stars and to study their distribution on the Milky Way bulge. 

Methods. Eor each star in the tile with more than 25 epoch (~90% of the objects down to ~ 17.0 mag), the standard deviation and 
y- test were calculated to identify variable candidates. Periods were determined using the analysis of variance. Objects with periods 
in the RR Lyrae range of 0.2 < P < 1.2 days were selected as candidate RR Lyrae. They were individually examined to exclude false 
positives. 

Results. A total of 1.5 sq deg were analyzed, and we found 39 RR Lyr stars, 27 of which belong to the ab-type and 12 to the c-type. Our 
analysis recovers all the previously identified RR Lyrae variables in the field and discovers 29 new RR Lyr stars. The reddening and 
extinction toward all the RRab stars in this tile were derived, and distance estimations were obtained through the period-luminosity 
relation. Despite the limited amount of RR Lyrae stars studied, our results are consistent with a spheroidal or central distribution 
around ~ 8.1 and ~ 8.5 kpc. for either the Cardelli or Nishiyama extinction law. Our analysis does not reveal a stream-like structure. 
Nevertheless, a larger area must be analyzed to definitively rule out streams. 

Key words. Galaxy: bulge - Galaxy: stellar content - Variable stars: RR Lyrae - surveys - Infrared: stars - Stars: horizontal-branch 


1. Introduction 

The structure and evolution of the Galactic bulge is far from 
completely explored. The high complexity of the Milky Way 
formation and the fact that we are embedded in the disk of 
the Galaxy, impede direct observations because of the extensive 
layer of dust and gas of the Galaxy in this line of sight. In the era 
of large-surveys, much effort have been made to reveal its struc¬ 
ture, taking advantage of a wide variety of near-, mid- and far-IR 
filters to bypass the layer of dust that does not permit deeper ob¬ 
servations in the optical bands (e.g., Freudenreich| 1 996} Skrut- 
skie et al.|2006[ Wright et al.|2010| l. 

In this context, the VISTA Variables in the Via Lactea (VVV) 
Survey aims to probe the internal structure of the Milky Way 
bulge through the near-IR (ZYJHK^) filters, using as its main 
distance indicators the red-clump giants and pulsating variable 
stars (Minniti et al. 2010| l, such as RR Lyrae stars (hereafter 
RR Lyr), classical Cepheids, anomalous Cepheids and Miras, 
and semiregular variables. In particular, RR Lyr stars have been 


19951 IGratton & Carrettal 

T99^ IMcWilliam & Zoccali |2010l 

Soszynski et al. |2011[ De 

cany et al. 201311. RR Lyr stars are 


old and low-mass stars O.TMq, |Smith||2b04l l that are in the 
horizontal branch (HB) stage (helium-burning phase) and have 
a very tight period-luminosity {P - L) relation in the near-IR 
bands (|Longmore et al. 1990[ Catelan et al. |2004|l. This rela¬ 


tion is a powerful tool that has previously provided distances to 


Galactic globular clusters (|Longmore et a .||1990||Del Principe 

et al. 

2005 2006[ Sollima et al.|20061 Coppola et al.|201 IJl, ex- 

traga 

actic globular clusters, dSph galaxies (jDalTOra et al. 2004} 

Pietrzynski et al.|2008l|Borissova et al.|2009|l, and the Galactic 

center (|Carney et al. 1995 

Groenewegen et al. 20081 Dekany 

|et al.|20131l. Dekany et al. ( 

2013|l searched for near-IR counter- 


* Based on observations taken within the ESO VISTA Public Survey 
VVV, Programme ID 179.B-2002. 


parts of the RR Lyr stars that were found by the OGLE III survey 
( |Soszynski et al. 20111 in the Galactic bulge in an area limited 
by 2° < \b\ < 7° and \€\ < 10°. They derived the distance to 
the Galactic center and discovered that the RR Lyr trace a more 
spheroidal shape than the red-clump stars, which suggests that 
the Milky Way has a composite bulge. 
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In this paper, we have selected a low density field to con¬ 
duct a variability search. We find new RR Lyr stars in the VVV 
tile b201, which allows us to measure precise reddening and dis¬ 
tances to the bulge in this direction. The tile b201 samples the 
spatial distribution of RR Lyr in the older outer bulge, b > 4° 


away from the Galactic center region analyzed by Dekany et al. 
( |2013| l. Our motivation also is to update the known RR Lyr stars 
present at these longitudes {{ ~ -10°) that are covered by the 
VVV Survey area. 


2. Observations 

The VVV Survey is being carried out with the 4m Visible and 
Infrared Survey Telescope for Astronomy (VISTA) located at 
ESO’s Cerro Paranal Observatory in Chile, using the wide-held 
VISTA InfraRed Camer a (VIRCAM; [Dalton et al.|2006l|Emer-| 
Ison & Sutherland|2010|l. Photometric catalogs are based on the 


VISTA system, for which 2MASS ( jCutri et'aL 


20031 coordi¬ 


nates are automatically produced by the Cambridge Astronom¬ 
ical Survey Unit (CASU|^ They are publicly available through 
the VISTA Science Archive (VSA^ The pixel scale and held 
of view (EoV) are 0.34" pix“^ and 1.64 deg^. This EoV repre¬ 
sents one “tile”; the entire VVV observations comprise of 348 
tiles (196 covering the bulge and 152 and the southern Galac¬ 
tic plane). The VVV observation schedule includes single-epoch 
near-IR photometry in ZYJHK^ bands and a variability cam¬ 
paign in with up to a total of 100 epochs planned by the end 
of the scheduled observing time (1929 hours, [Minniti et al.|2010[ 
jCatelan et al. [20116) 1. Technical details of the telescope and ob¬ 
servation strategy can be found injMinniti et al.|(|2010|l andjSaito 
|etar] ( |20T2al l. 

In this study, we used observations of tile b201 ({,b ^ 
-9°,-9°; Pig. [T]). This tile is the bulge held most distant from 
the Galactic center observed by the VVV, and it has the lowest 
stellar density, as shown in Pig. 2 of jSaito et al. ( 2012b| ). The 
VVV has observed this tile 36 times between April 2010 and 
July 2013. A total of 394047 point sources with Ks < 17.0 mag 
were detected, of which 353 935 have more than 25 epochs of 
observations, which is the number we require to include a source 
in this study. 

The main reason to choose this tile was its stellar density, 
lower than other observed VVV bulge tiles. Specihcally, for a 
cutoff of 25 points per light curve, there were 353 935 time se¬ 
ries to be analyzed in this tile of the total 394047 (~ 10 % ex¬ 
cluded). This number is similar to that of the adjacent tile b202 
with 408 858 sources in band. Another reason for analying 
this tile is its low extinction with values of <0.1 mag in cer¬ 
tain lines of sight, or equivalently. Ay < 0.35 mag (Catelan et al 


201 lb[ Gonzalez et al.|20rT] 2012|l, w hich allowed the previous 
detection of RR Lyr stars by Swope|([l942|l,|Po^nsen|(|1954|l and 


Kooreman ( 1966|l. Even though it is not completely free of se- 
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Fig. 1. VVV Survey bulge region, shown in Galactic coordinates with 
its individual numbered tiles. Tile b201 is highlighted in the bottom 
right comer. The Galactic center is contained in tile b333. Adapted from 
[Catelan et al.|j2011b| >. 


(crKs) test for each time series. The cut-off for consider¬ 

ing a non-variable star was determined by simulating a random 
noise distribution, as described by Carpenter et al. ( 200l| ; 
this resulted in - 2. This constraint places a lower limit of 
Aa:, ~ 0.03 mag at /Tj ^ 13 mag on the amplitudes of the 
RR Lyr to be detected. We estimate, based in the w Cen RRc stars 
Kj-amplitude distribution (Navarrete et al. 2014, subm.), that at 
Ks ~ 15.2 mag our search procedure has a variability detection 
efficiency of ~9Q%, which drops to ~50% at Ks ~ 16.5 mag. 
These candidates were analyzed by the analysis of variance 
(AoV) statistis ( |Schwarzenberg-Czerny| 1989) 1 to determine their 
periods. Periods were restricted to the range in which RR Lyr are 
typically found: 0.2 < P(days) < 1.2. After this step, the light 
curves of 167 candidates were visually inspected, and 42 objects 
were foundto be as RR Lyr. Of the remaining 125 light curves, 
32 were clearly variables but without a clear variability class, 
and we classified the remaining 83 as eclipsing-binary candi¬ 
dates. Better phase coverage is needed to classify these variables. 
Pinally, ten variable-star candidates as selected by the algorithm 
were determined to be non-variables that showed period-aliasing 
of ail day and were therefore excluded from further analysis. 
Pig. [^presents the magnitude-variability diagram of all sources 
with more than 25 epochs of observation, in which the RR Lyr 
stars are identified. 


lection biases, the low density implies reduced crowding and a 
more accurate and deeper photometry than for the other fields in 
the VVV survey area. Purthermore, the lower extinction in this 
area minimizes the reddening uncertainties. 

2.1. Identifying and classifying variable sources 

A semi-automated classification scheme was developed to se¬ 
lect the variable stars. Pirst we computed the standard deviation 


http://casu.ast.cam.ac.uk/vistasp/ 
http://horus.roe.ac.uk/vsa/ 


3. Results 

3.1. Catalogue of RR Lyr stars 

To characterize our sample of RR Lyr stars. Table [T] contains 
the number of each variable type found and the mean, longest 
and shortest periods of each population. These values define the 
Oosterhoff and Bailey type of the RR Lyr. Additionally, our 
RR Lyr candidates lie in the magnitude range between 13 < 
Ks (mag) < 16.5, resulting in a distance between 5.2 < d (kpc) < 
13.6^ _ _ 

[Bailey I ( 1902| ) classified the RR Lyr (in that epoch called 
“cluster variables”) in three groups, the a-, b- and c-type by its 
period and light-curve shape (amplitude of variation). The first 
two groups were merged later into the ab-type known today. This 
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Fig. 2. Magnitude-variability diagram for the sources detected in b201 
with more than 25 epochs (levels of blue). The contribution of saturated 
(Ks < 12.0 mag) and faint (K^ > 16.5 mag) stars are clearly visible 
in the crKs value. RRab are shown as red stars and RRc stars as black 
circles. 
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Fig. 4. Hess diagram of the tile b201. The disk MS and bulge RGB 
are easily recognizable from this field and similar to the ones derived 
in|Saito et al.|l|2012b|l. The symbols of the RR Lyr are the same as in 
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instability strip locus can be inferred over the disk MS in a color 
range between 0.1 < (J - < 0.6, according to the position of 

the RR Lyr in the diagram. 

The period distribution of the RRab and RRc stars is divided 
by a gap in P ~ 0.4 days, as shown in the b ottom panel of Fig. 
which reflects the distribution derived by |Kunder & Chabo}^ 
P009| l in V and R bands. These authors analyzed the Ooster- 
hoff type of RR Lyr in the Galactic bulge. This near-IR diagram 
shows a different slope for the RRab in the optical wavelengths. 
It is more horizontal, similar to that reported by |Gavrilchenko| 
et al. (|2014 1 and Navarrete et al. (2014, subm.) using mid- and 


near-IR data. 


Fig. 3. Top: Bailey diagram for the 39 RR Lyr variable stars. As ex¬ 
pected, there is a clear separation between the amplitude and period 
among RRab (upper right) and RRc stars (lower left), following the 
same symbols as in Fig.|^ Bottom: Period histogram of the RR Lyr. 

Table 1. Summary of the periods (in days) of the RR Lyr found in b201. 


Variability type 

Quantity 

(P) 

(days) 

P . 

^ mm 

(days) 

P 

-* max 

(days) 

RRab 

27 

0.578 

0.448 

0.787 

RRc 

12 

0.312 

0.202 

0.379 


relation of period and amplitude (today called the Bailey dia¬ 
gram) is shown in the upper panel of Fig. 

J and H photometry taken during 2010 and the fTs-mean 
magnitudes are shown in Table A brief discussion of the er¬ 
ror introduced in the RR Lyr color {J - K^) and (H - is 
presented in Sect. |3.5| The RR Lyr stars are shown on a Hess 
color-magnitude diagram of the whole tile in Fig. The bulge 
red-giant branch (RGB) and the disk main-sequence (MS) stars 
are presented in the Hess diagram, with the label aligned in the 
direction of each branch. The superposition of a young popu¬ 
lation (the disk) hides the MS turn-off and part of the Galactic 
bulge HB, where most of our RR Lyr belong. However, the bulge 


3.2. Previously known RR Lyr stars 


As a test of consistency, we compared our RR Lyr positions and 
periods with those found in the literature. There are already ten 
variables in the International Variable Star Inde)([^(VSX) catalog 
classified as RR Lyr in the analyzed area. 

We recovered all ten previously known variables and cor¬ 
rectly classified them as RR Lyr. These ten RR Lyr are listed 
in Table 3] along with their General Catalogue of Variable Stars 
(GCVS; Samus et al.|2009| l designations, periods, and astromet¬ 
ric offsets. In three cases (MO CrA, V397 CrA and V467 CrA), 
periods different from those in the literature were identified. MO 
CrA was classified, without an assigned period, as a “cluster 
variable” by Swope|| ( l942|l. V 397 CrA and V467 CrA were pre¬ 
viously observed by Ponsen ( |1954| l and Kooreman ( 1966| l, re¬ 
spectively. |Ponsen ]( |1954| l classified V397 CrA as an RRc, but 
with only 16 points in the light curve. Kooreman ( 1966| l classi¬ 
fied V467 CrA as an RRab, remarking that this star was very 
faint, with a large dispersion in magnitude in the 120 avail¬ 
able plates. V475 CrA and V493 CrA do not have counterparts 
within our initial 2"matching radius, but RR Lyr with near iden¬ 
tical periods were identified at 17"and 31", respectively, from 
their VSX positions. According to CASH, VVV astrometry is 


^ http://www.aavso.org/vsx 
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Table 2. Color information of RRab and RRc stars found in tile b201. We list the VVV ID, right ascension (J2000), declination (J2000), magnitude- 
weighted mean from the Fourier fit and (J - K^) and (H - K^) colors from the first Kj-band epoch. 


VVV ID (RRab) 

RA(J2000) 

DEC(J2000) 



(H-KX 

VVV J2712426.16-421241.6 

18:05:37.7 

-42:12:41.6 

14.59 

0.54 

0.22 

VVV J2712711.57-421314.0 

18:05:48.8 

-42:13:14.0 

14.23 

0.65 

0.29 

VVV J2712912.05-420223.7 

18:05:56.8 

-42:02:23.7 

14.23 

0.44 

0.17 

VVV J2711512.06-415330.1 

18:05:00.8 

-41:53:30.1 

14.00 

0.51 

0.22 

VVV J2703434.23-413322.1 

18:02:18.3 

-41:33:22.1 

13.35 

0.57 

0.24 

VVV J2703536.01-412829.4 

18:02:22.4 

-41:28:29.4 

13.97 

0.42 

0.20 

VVV J2703343.00-412731.0 

18:02:14.9 

-41:27:31.0 

14.45 

0.60 

0.25 

VVV J2711523.24-413619.8 

18:05:01.5 

-41:36:19.8 

13.78 

0.64 

0.27 

VVV J2710804.03-41324L7 

18:04:32.3 

-41:32:41.7 

14.51 

0.50 

0.21 

VVV J2704944.51-412248.6 

18:03:19.0 

-41:22:48.6 

13.36 

0.34 

0.16 

VVV J2714215.15-414058.6 

18:06:49.0 

-41:40:58.6 

14.29 

0.64 

0.25 

VVV J2710759.32-411700.5 

18:04:31.9 

-41:17:00.5 

13.39 

0.41 

0.20 

VVV J2714823.59-41312L7 

18:07:13.6 

-41:31:21.7 

14.44 

0.26 

0.15 

VVV J2711644.06-411744.6 

18:05:06.9 

-41:17:44.6 

13.84 

0.50 

0.21 

VVV J2710158.99-410727.3 

18:04:07.9 

-41:07:27.3 

14.22 

0.42 

0.20 

VVV J2704656.26-421805.3 

18:03:07.7 

-42:18:05.3 

15.22 

0.64 

0.26 

VVV J2701149.18-420358.3 

18:00:47.3 

-42:03:58.3 

13.22 

0.59 

0.25 

VVV J270510L04-405300.9 

18:03:24.1 

-40:53:00.9 

13.58 

0.44 

0.18 

VVV J2711836.69-405946.8 

18:05:14.4 

-40:59:46.8 

14.24 

0.56 

0.25 

VVV J2705334.00-421303.7 

18:03:34.3 

-42:13:03.7 

15.04 

0.67 

0.28 

VVV J2702317.06-420003.9 

18:01:33.1 

-42:00:03.9 

14.14 

0.50 

0.23 

VVV J2705622.48-420853.7 

18:03:45.5 

-42:08:53.7 

14.08 

0.66 

0.29 

VVV J2711018.90-421253.4 

18:04:41.2 

-42:12:53.4 

14.03 

0.60 

0.27 

VVV J2712342.04-421647.8 

18:05:34.8 

-42:16:47.8 

14.13 

0.28 

0.15 

VVV J2712855.63-421832.3 

18:05:55.7 

-42:18:32.3 

14.46 

0.40 


VVV J2710638.97-420755.7 

18:04:26.6 

-42:07:55.7 

14.59 

0.43 

0.19 

VVV J2702400.61-415026.5 

18:01:36.0 

-41:50:26.5 

14.72 

0.34 

0.18 

VVV ID (RRc) 

RA(J2000) 

DEC(J2000) 


(J-K,) 

(H-K,) 

VVV J2711026.32-415037.7 

18:04:41.7 

-41:50:37.7 

15.39 

0.29 

0.13 

VVV J2702945.91-413239.9 

18:01:59.0 

-41:32:39.9 

12.73 

0.33 

0.09 

VVV J2710200.43-413446.0 

18:04:08.0 

-41:34:46.0 

12.84 

0.05 

0.02 

VVV J2705846.71-411723.2 

18:03:55.1 

-41:17:23.2 

14.63 

0.43 

0.17 

VVV J2705404.77-422910.0 

18:03:36.3 

-42:29:10.0 

14.36 

0.13 

0.08 

VVV J2713458.02-412608.5 

18:06:19.9 

-41:26:08.5 

13.63 

0.10 

0.06 

VVV J2711142.24-411611.0 

18:04:46.8 

-41:16:11.0 

14.79 

0.18 

0.10 

VVV J2714054.14-412637.6 

18:06:43.6 

-41:26:37.6 

14.33 

0.01 

0.02 

VVV J2711803.26-411635.8 

18:05:12.2 

-41:16:35.8 

11.70 

0.05 

0.01 

VVV J2710857.12-410712.8 

18:04:35.8 

-41:07:12.8 

14.90 

0.37 

0.17 

VVV J2702329.86-415953.8 

18:01:34.0 

-41:59:53.8 

16.51 

0.61 

0.22 

VVV J2704051.15-420503.6 

18:02:43.4 

-42:05:03.6 

13.07 

0.12 

0.10 


Notes. Typical (J - H) and (H - K^) photometrical errors are 0.05 mag. //-magnitude not available. 


accurate to ^50 mas and the similarity in periods make it very 
likely that V475 CrA and V493 CrA are counterparts to VVV 
J2705622.48-420853.7 and VVV J2714823.59-413121.7. 


3.3. RR Lyr star with the shortest period in our cataiog 

We found that the RRc candidate VVV J2705846.71-411723.2 
has a very short period of P = 0.202 days. In contrast, the 
RRc with the shortest period found in the OGLE III sample has 
P = 0.237 days ( |Soszyhski et al.|2011| l. This variable warrants 
a more detailed study to confirm that it is an RRc instead of an 


eclipsing binary or a long-period SX Phe, as reported by Cohen 
|& Sarajedini| ( |2012] l. 


3.4. Fourier coefficients 

The complete catalog of light curves is shown in Figs. |^and|^ 
for the RRab and RRc types. Following an approach similar to 
that of |Dekany et al.| ( |2013| l, the light curves of the RRab and 
RRc stars were fitted with Fourier series compound of sines, up 
to sixth order using the direct Fourier fitting (DFF) method of 
|Kovacs & Kupi|p007| l. This order was the highest possible cal¬ 
culated, but in some cases the returned order was lower to pre¬ 
vent over-fitting of the light-curve. This method allows us to re¬ 
cover the mean /fs-magnitude plus the Fourier coefficients (A 21 , 
(f> 2 i, ..., Agi, which are shown in Table|^ 

The photometric iron abundance of RR Lyr stars is typically 


found by analyzyng of V- and /-band light-curves (e.g., Jurcsik 


|& Kovac^|1996[ |Smolec|[2005|). However, since there is no re- 
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Table 3. Results of the match with the VVV tile b201 data and the VSX catalog. 


VVV ID 

GCVS 

Fvvv 

(days) 

Frit. 

(days) 

d 

n 

VVV J2705101.04-405300.9 

MO CrA 

0.657005 


0.475 

VVV J2712855.63-421832.3 

V397 CrA 

0.491175 

0.3293815 

0.771 

VVV J2701149.18-420358.3 

V463 CrA 

0.604052 

0.6040585 

0.139 

VVV J2702329.86-415953.8 

V467 CrA 

0.359989 

0.4480160 

1.931 

VVV J2703434.23-413322.1 

V469 CrA 

0.678028 

0.6780236 

0.528 

VVV J2705622.48-420853.7 

V475 CrA 

0.511933 

0.5119430 

17.41 

VVV J2710759.32-411700.5 

V482 CrA 

0.541709 

0.5417140 

2.725 

VVV J2712342.04-421647.8 

V483 CrA 

0.485059 

0.4850490 

0.105 

VVV J271271L57-421314.0 

V486 CrA 

0.460161 

0.4601694 

0.573 

VVV J2714823.59-41312L7 

V493 CrA 

0.519262 

0.5194291 

31.26 


lationship between the Fourier parameters and iron content al¬ 
ready established in the near-lR, it cannot be determined individ¬ 
ually for our variables. Hence, we used the value of|Pietrukowicz 


et al. (2012|l, who found that RR Lyr trace the distribution of the 


metallicity in the Galactic bulge with a sharply peaked distribu¬ 
tion, centered on -1.02 dex with a dispersion of 0.25 dex. This 
metallicity complete agree with the spectroscopic value found 
by Walker & Terndru^ ( 1991| l of {{FelH}) - -1.0 dex with a 
dispersion of 0.16 dex. They analyzed 59 RR Lyr in Baade’s 
window. Taking a fixed metallicity does not significantly affect 
the distance estimation; it imparts a mean of 0.20 kpc (~ 3% 
at 8 kpc) uncertainty to distances derived using the P - L rela¬ 
tion, taking into account the extreme values on the dispersion of 
[Fe/H]. 


2014| l is collecting. The next step was to calculate the intrinsic 
color of the sample, that is (J - Kfi)o, analogous to the difference 
in absolute magnitude of the two filters. This modifies the latter 
reddening equation into 


E{J - K,) = (7 - K,) -{J- K,)o = (7 - K,) - (Mj - Mr,) (4) 

for 7 and bands. The absolute magnitudes were calculated 
using the P - L relation given in|Alonso-Garcfa et al.|(|2014)l. 


Ms = -0.6365 - 2.347 log F -I- 0.1747 log Z, (5) 


3.5. RR Lyr stars: reddening estimation 

Generally, the reddening can be expressed as the difference of 
the measured and intrinsic color between two filters. 


E(J - K,) ^ (J - K,) - (J - K,)o, (1) 

in this case for 7 and The first term of the equation is 
limited by the phase coverage and the number of points of each 
variable in multiple wavelengths. The VVV, as a variability sur¬ 
vey, provides time-series exclusively in band. For 7 and El 
band, only one observation has been taken during the first year 
of operation, which limits the precision of (J - K^) and (H- K^). 
The colors (J-H) and were calculated using the faintest 

point of the light curve according to the Fourier fit, and taking 
into account the ratio of amplitudes between 7, H and de¬ 
scribed in |Feast et al.| ( |2008| l, as follows: 

AH =-1-0.11-H 1.65(A/:s-0.18), (2) 


A7 =-0.02 h-3.60(AHs-0.18). (3) 

This approximation has a larger error for the shorter-period 
RRab (higher amplitudes) with a maximum of +0.2 mag and 
+0.1 mag for the 7 and H band, respectively. However, for Z 
and Y bands (also available in the VVV Survey), we are unable 
to conduct the same analysis because we lack detailed near-lR 
variability studies. This problem can be solved in the future with 
a database of well-defined near-lR template light-curves that the 
VVV Templates Project ( |Catelan et al.|[2011a[ [Angeloni et al.| 


Mk, = -0.6365 - 2.347 logP + 0.1747 logZ, 


( 6 ) 


with cTj = 0.17 mag and = 0.13 mag the systematic er¬ 
rors in 7 and bands, respectively. The errors were calculated 
as the difference between the average magnitudes, published by 
[Del Principe et ah] ( |2005[ ), and the theoretical values obtained 
using the P - L relations listed above. The theoretical propa¬ 
gated error of (7 - Hs)o is ~ 0.21 mag, which introduces a dif¬ 
ference of 0.14 mag in the distance moduli. These relations are 
crucial for determining the distances to our RR Lyr stars, be¬ 
cause it is a very accurate solution to calculate the reddening 
corrected colors, especially in the near-lR ([Catelan et al.|2004 


We used the same approach as Pietrukowicz et al. ( |2012[ ) to ca 
culate the logZ = [Fe/H] - 1.765, based on a solar metallicity 
of Zq = 0.017 ( [Catelan et al.||20D4| ). This relation was estab¬ 
lished for fundamental-mode pulsators, therefore we must ex¬ 
clude the RRc stars of future analyses. Finally, these colors give 
us an estimate of the individual reddening of the RRab stars, 
which is shown in Tabl e [5] We compa red these reddening values 
with those obtained by jKunder et al.| ( f2008| l in the same Galac¬ 
tic latitude as our values ({,b ~ -5°,-9°) from the Galactic 
bulge fields of the MACHO survey. The mean converted value 
of E{J - Ks) at that latitude is ~ 0.20 mag, with a dispersion of 
0.04 mag, which is consistent with our measurements within the 
errors. In some cases, our estimates are slightly higher than those 
of jKunder et al.|p008| l, which may be due to larger distances or 
larger errors in the determination of the color (7-Hs) in our light 
curves. 


3.6. Distance to the RR Lyrae stars 

We calculated the individual distances to the RRab stars using 
the P - L relations for the band. This distance can be ex- 
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Fig. 5. 30 RRab stars found in tile b201, sorted by increasing period. The solid line represents the Fourier decomposition of each RR Lyr using the 
DFF method. The internal name and period are indicated at the top of each light-curve. Error bars are plotted, but are generally smaller than the 
point sizes. 
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VWJ2705846.71-411723.2 (P = 0.202) 



VWJ2714054.14-412637.6 (P = 0.251) 



VWJ2711026.32-415037.7 (P = 0.252) 



VWJ2710857.12-410712.8 (P = 0.285) 





VWJ2713458.02-412608.5 (P = 0.331) 



VWJ2702945.91-413239.9 (P = 0.333) 



VWJ2711142.24-411611.0 (P = 0.359) 



VWJ2702329.86-415953.8 (P = 0.36) 



VWJ2711803.26-411635.8 (P = 0.371) 



VWJ2710200.43-413446.0 (P = 0.379) 



Fig. 6. Same as Fig.|^ but for the RRc stars. 


pressed in terms of the absolute and the extinction-corrected K^- 
band magnitudes. 


logf? = 1 + 0.2(/:,,o - Mk,), 


(7) 


with R the individual distance in pc to our RRab stars and 
converting the color excess into extinction using the law of 
ICardelli et ari ( |T989] l, 


Aks = 0.73E(J - K,). 


( 8 ) 


Figure 1^ shows the histogram of the distances calculated for 
the individual RR Lyr. A strong agreement with previous results 
is recovered even at this high latitude ~ -10°) with more 
than the 50% of the variables in the central kiloparsec accord¬ 
ing to a Gaussian fit with r/mean = 8.14 kpc and a dispersion of 



(Td = 1 .56 kpc. Adopting the reddening law of Nishiyama et al. 
( |2009| l does not significantly change the distance distribution of 
our RRab stars; it adds 0.4 kpc to the mean, but keeps the dis¬ 
persion unchanged. 

The results of this analysis, including reddening, extinction, 
and distances with their error are shown in Table |5] The uncer¬ 
tainties consider the error propagation on the color, reddening, 
extinction, and distance of each individual RRab star. Another 
visualization of the derived distances, assuming a constant i and 
b coordinate, is shown in Fig as a cone view. This approach 
is useful in identifying groups or streams of variable stars in our 
FoV, which is not possible with their spatial distribution {£, b) 
alone, as also shown in Fig. Although there may be a group 
of RR ah stars at ~ 7 kpc, these stars do not cluster in any 
other parameters, for example, in the period-amplitude diagram 
or spatially, and hence we are hesitant to consider them as a 
group. 


Fig. 7. Distance to the sample of RR Lyr stars in tile b201 in a 1.2 kpc 
binning size. The solid (D13) and dashed lines (PI 1) represent the val¬ 
ues of the Galactic center distance given in |Dekany et al.|j2013) and 
[Pietrukowicz et al!]P012| >. 


An interesting topic of the RR Lyr studied in large areas is 
the possibility that they can trace an old stream-like structure 
( |Sesar et al.|2013[) in a ddition to the spheroidal shape described 
by |Dekany et al. ( 2013| l. Unfortunately, the small number of vari¬ 
ables in this single tile prevent us from making any strong state¬ 
ments about the presence of streams in this direction. 
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Table 4. Period and Fourier coefficients for the RR Lyr in tile b201. 


VVV ID (RRab) 

P (days) 

^1 

1^1 

A 21 

1^21 

A 31 

^31 

A 41 

<^41 

A 51 

1^51 

Aei 

<I>61 

VVV J2702317.06-420003.9 

0.448043 

0.108 

0.146 

0.442 

3.361 

0.156 

0.038 

- 

- 

- 

- 

- 

- 

VVV J2712711.57-421314.0 

0.460161 

0.127 

0.595 

0.454 

3.768 

0.336 

1.103 

0.179 

-1.599 

0.156 

1.955 

0.132 

-0.353 

VVV J27I2342.04-421647.8 

0.485059 

0.098 

1.042 

0.336 

3.942 

0.109 

1.388 

- 

- 

- 

- 

- 

- 

VVV J2712855.63-421832.3 

0.491175 

0.122 

5.804 

0.376 

-8.542 

0.196 

-17.03 

- 

- 

- 

- 

- 

- 

VVV J2712426.16-421241.6 

0.510823 

0.120 

3.069 

0.403 

-2.386 

0.228 

-4.592 

- 

- 

- 

- 

- 

- 

VVV J2705622.48-420853.7 

0.511933 

0.137 

5.191 

0.413 

-8.957 

0.174 

-11.045 

0.123 

-20.52 

0.117 

-23.676 

0.099 

-21.225 

VVV J27I0158.99-410727.3 

0.513884 

0.128 

2.769 

0.295 

-2.487 

0.157 

-4.873 

0.116 

-7.579 

- 

-13.846 

- 

- 

VVV J2714215.15-414058.6 

0.516757 

0.116 

1.954 

0.393 

-2.456 

0.183 

-5.302 

0.098 

-1.648 

- 

-9.772 

- 

- 

VVV J2714823.59-413121.7 

0.519262 

0.105 

2.818 

0.317 

-2.442 

0.118 

-4.912 

- 

- 

- 

- 

- 

- 

VVV J2705334.00-421303.7 

0.523466 

0.118 

0.447 

0.432 

4.396 

0.310 

2.112 

0.175 

0.423 

0.197 

3.431 

- 

- 

VVV J2703343.00-412731.0 

0.535809 

0.117 

5.564 

0.321 

-8.513 

0.258 

-11.026 

0.102 

-19.275 

0.072 

-27.465 

- 

- 

VVV J2710804.03-413241.7 

0.536988 

0.114 

4.762 

0.355 

-8.314 

0.159 

-10.869 

0.141 

-13.575 

0.131 

-21.91 

0.121 

-24.132 

VVV J2710759.32-411700.5 

0.541709 

0.113 

1.205 

0.378 

-2.069 

0.254 

2.238 

0.111 

0.374 

- 

- 

- 

- 

VVV J2711836.69-405946.8 

0.543833 

0.126 

2.396 

0.290 

-1.936 

0.173 

-4.095 

0.187 

-7.041 

0.125 

-9.401 

- 

- 

VVV J2711523.24-413619.8 

0.548459 

0.113 

4.514 

0.407 

-8.315 

0.185 

-10.374 

0.177 

-12.632 

0.123 

-20.615 

- 

- 

VVV J2710638.97-420755.7 

0.560748 

0.107 

2.260 

0.330 

-2.067 

0.159 

-4.570 

0.114 

-7.347 

- 

- 

- 

- 

VVV J2702400.61-415026.5 

0.560889 

0.108 

2.005 

0.325 

-2.066 

0.239 

-4.529 

0.131 

-6.443 

0.069 

-8.248 

- 

- 

VVV J2701149.18-420358.3 

0.604052 

0.122 

3.905 

0.387 

-2.170 

0.255 

-11.105 

0.235 

-13.759 

0.131 

-16.498 

- 

- 

VVV J2711512.06-415330.1 

0.625496 

0.122 

2.897 

0.375 

-2.041 

0.271 

-4.298 

0.130 

-6.362 

0.055 

-14.464 

0.055 

-15.268 

VVV J2711018.90-421253.4 

0.638694 

0.136 

1.014 

0.404 

4.155 

0.304 

1.629 

0.218 

-1.227 

0.172 

-3.558 

0.089 

-0.497 

VVV J2712912.05-420223.7 

0.642877 

0.109 

4.470 

0.266 

-8.214 

0.130 

-9.892 

0.091 

-17.681 

- 

- 

- 

- 

VVV J2705101.04-405300.9 

0.657015 

0.116 

0.131 

0.238 

4.176 

- 

- 

- 

- 

- 

- 

- 

- 

VVV J2704656.26-421805.3 

0.660970 

0.133 

4.211 

0.366 

-8.292 

0.152 

-10.250 

0.137 

-13.372 

0.144 

-20.435 

- 

- 

VVV J2703434.23-413322.1 

0.678028 

0.121 

3.676 

0.386 

-2.160 

0.234 

-10.655 

0.114 

-12.769 

0.074 

-14.812 

- 

- 

VVV J2704944.51-412248.6 

0.746910 

0.113 

5.754 

0.241 

-7.914 

0.126 

-16.045 

0.102 

-18.418 

0.104 

-28.086 

- 

- 

VVV J2703536.01-412829.4 

0.756627 

0.126 

2.947 

0.372 

-2.210 

0.214 

-4.254 

0.124 

-6.989 

0.063 

-8.932 

- 

- 

VVV J2711644.06-411744.6 

0.786726 

0.122 

4.691 

0.207 

-8.180 

0.119 

-10.891 

0.114 

-13.366 

0.110 

-21.196 

- 

- 


VVV ID (RRc) 

P (days) 


<^i 


4>21 

A 31 

<p3l 

VVV J2705846.71-411723.2 

0.201895 

0.089 

6.168 

0.109 

-7.217 

0.101 

-14.489 

VVV J2714054.14-412637.6 

0.251285 

0.033 

2.613 

- 

- 

- 

- 

VVV J2711026.32-415037.7 

0.252489 

0.113 

5.647 

- 

- 

- 

- 

VVV J2710857.12-410712.8 

0.285146 

0.124 

4.937 

0.066 

-7.98 

0.057 

-13.387 

VVV J2704051.15-420503.6 

0.302654 

0.076 

0.275 

- 

- 

- 

- 

VVV J2705404.77-422910.0 

0.322360 

0.055 

5.595 

- 

- 

- 

- 

VVV J2713458.02-412608.5 

0.330849 

0.060 

0.422 

- 

- 

- 

- 

VVV J2702945.91-413239.9 

0.332772 

0.019 

5.313 

0.236 

-4.388 

- 

- 

VVV J2711142.24-411611.0 

0.358777 

0.080 

3.571 

0.255 

-2.865 

- 

- 

VVV J2702329.86-415953.8 

0.359989 

0.105 

0.209 

0.409 

5.568 

- 

- 

VVV J2711803.26-411635.8 

0.370611 

0.044 

1.560 

- 

- 

- 

- 

VVV J2710200.43-413446.0 

0.378936 

0.019 

5.662 

0.220 

-9.100 

- 

- 


4. Summary 

A search for RR Lyr stars was performed in the VVV tile b201, 
one of the outermost bulge regions covered by the VVV Sur¬ 
vey. A total of 39 RR Lyr were found, of which 27 and 12 be¬ 
long to the ah- and c-type, respectively. For the total sample, 
coordinates, periods, amplitudes, and near-lR colors were pre¬ 
sented. Using the International Variable Star Index, we found 
counterparts of the ten previously known RR Lyr in the field, 
in agreement with past measurements. We found differences in 
only three of them, for which we suggested a revised period. 
From the analysis of the RRab sample, individual reddening, ex¬ 
tinction, and distance were obtained, based on line-of-sight ex¬ 
tinction values given by Cardelli et al.| ( |1989[ ) and the R - L re¬ 
lation of |Alonso-Garcfa et al.| ( |2014| l. The distribution of these 
variables on the Milky Way follows a central distribution around 
~ 8.1 and ~ 8.5 kpc. using the Cardelli or Nishiyama extinc¬ 
tion law, respectively. This value is consistent with the results 
of Dekany et al. ( 2013| l, placing the center of the distribution at 

~8.3 kpc. An analysis of a larger area around f -10° will be 

conducted to complement the studies of the inner Galactic bulge 
that are already published. 
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Table 5. Reddening, extinction and distances computed for the RR Lyr found in tile b201. 


VVV ID (RRab) 

E{J-K,) 

(mag) 

(mag) 

Distance 

(kpc) 

VVV J2712426.16-421241.6 

0.33 

± 

0.13 

0.23 
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0.31 
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VVV J2703536.01-412829.4 
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± 

0.09 

8.79 
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0.28 

VVV J2703343.00-412731.0 

0.38 
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0.22 
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Fig. 8 . Left: Distance to the sample of RR Lyr stars, between 5 and 10 kpc, in a cone view. We did not find any apparent cluster or group in the 
RRab due to streams, except for the variables in the Galactic center. The solid and dashed lines have the same meaning as in Fig. 0 Right: Spatial 
distribution in Galactic coordinates of the 39 RR Lyr found in tile blOl. 
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